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SUMMARY 


Routine determination of inviscid subsonic flow fields about 
wing-body-tail configurations is state-of-the-art employing a Green’s 
function approach for numerical solution of the perturbation velocity 
potential eqmtion. This approach has been successfully extended 
into the hi^ subsonic subcritical flow regime and into the shock- 
free superaonic flow regime. However, it has not been successfully 
extended into the transonic flow regime. The present study develops 
a modified Green's function formulation, valid throughout a range of 
Mach numbers including transonic, that takes an explicit accounting 
of the intrinsic non-linearity in the parent governing partial 
differential equations. Some considerations pertinent to flow field 
predictions in the transonic flow regime are discussed. 
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IHTBODUCnON 


Potential flow aero^nasics has traditionally examined the solu- 
tion of the Laplacian of the incompressihle flow perturbation poten- 
tial function. This solution has found wide application in a variety 
of problems. Of particular interest here is its successful applica?- 
tion at subsonic and supersonic speeds by means of the Goethe3rt 
transformation. However, attempts to achieve a method of utilizing 
this solution at transonic speeds have proved unsuccessful. It is 
felt that this failure is due to shortcomings in the traditional 
formulation of the approximte potential equation at transonic speeds. 
The purpose of the present study is to evaluate the more comprehen- 
sive tensor form of the governing equation for irrotational potential 
flow, so as to extract a potentially useful Green's function solution 
form, valid throughout a Mach nunber range, specifically including 
transonic flow. 
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SYMBOLS 


function; coefficient matrix 
isentropic sound speed 
alternating tensor 
function; load vector 
function 

function; Green's function 
Green's tensor function 
unit vector triad 
differential operator on R 
scalar component of Mach tensor 
Mach number; tensor function 
normal vector 
pressure 

geodesic distance 


solution domain; tensor geodesic distance 
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perturbation velocity vector 
velocity vector 
Cartesian coordinate system 
rectangular Cartesian coordinate system 



ratio of specific heats 

Dirac delta function; Kronecker delta 

closxjre of R 

small parameter 

density 

singularity strength 
summation operator 
perturbation velocity potential 
total velocity potential 
function 


Subscripts 

» fi’eestream reference 


i.J,k 


t'^nsor indices 



Q 


0 


t 


normal to closure 
stagnation reference 
partial derivative by time 
gradient vector operator 
not a tensor index 


Superscripts 

sooirce point position vector 
-► vector 


unit vector 



REVIEW OF THE LINEARIZED PROBLEM FORMULATION 


CoBiputational solutions for three-dimensional, inviscid irrota- 
tional flows about practical aircraft configurations have achieved, 
a high state of capability. !Die theoretical fomulation is a Green's 
function solution of the Laplacian differential equation governing 
the distribution of velocity perturbation potential function. For 
incompressible flow, the continuity equation states 





( 1 ) 


Ihe second form in £q. (l) introduces the preferred tensor indicial 
notation to be eventxially required. The differential constraint of 
irrotationality is 


X U = 


®ijA,J 


= 0 


( 2 ) 


Equations (l) - (2) cwpletely describe the inconqjressible flow ve- 
locity vector. Furthermore, from Eq. (2), since the curl of the 
gradient vanishes identically, we can express U as 


U s s 

’k 


(3) 



Hence , the velocity field may be coogi^letely determined independent of 
the pressure t and neglecting momentum considerations (see reference 1, 
page 6), by solution of 



4 ,. 


kk 


(U) 


on the solution domain R, as determined by combining £q,s. (l) - (3)* 
£!q.. (U) is an elliptic boundary value description; the boundaxy condi- 
tions are typically Neumann, obtained by forming the inner product of 
Eq. (3) with a unit normal vector n on the closure, 3R, of R as 


4 • fi = U 

n 


(5) 


In Eq. (5), U is the local velocity efflux through 3R. Furthermore, 

tl 

from Eq. (3), the level of 4 is arbitrary and may be specified at a 
convenient value for some location on 3R. 

For small disturbance theory (see reference 2, chapter lU), the 
local flow velocity vector can be represented as 

U(x. ) ® U + u(x. ) (6) 

1 • 1 
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where u(x^ } ia the aoount hy which the local irrotstional velocity 
differs from the reference freeatrean value, U^. Assuming aligned 
parallel to the x axiSj and identifying the perturbation velocity po- 
tential. 




( 7 ) 


wo have 


u « a 


( 8 ) 


Since u describes an irrotational, inviscid flow, we obtain as before 


7^ d ■ 0 


( 9 ) 


Tiie Neumann boundary cwdition becomes, from Eq.. (8) 


-^(|) • n ■ (U - U,i) • n 


( 10 ) 


Ihe theory of Green's functions and potential theory are based 
v^on the divergence theorem 


/gV ■ fit . «do 


(n) 


If we let f - and respectively, insert each into Bq. (11) and 
siibtract, ve obtain Green's third identity. 


- 4>y^-4)]dT 



♦^3 • ndo 


( 12 ) 


Eq. (12) is the fo\mdation for a Green's function solution to Eqs. (9) - 
(10). Let tahe the specific form 


where 


^ 5 



(13) 





(ll*) 


We assume y^A - 0, and that A has no singularities in R. The function 

,-►,-1 

A can he used to adjust boundary values of g, and |r| is singular at 
X ■ X* , where x is the position vector of a field point and x' is the 
position vector to a source point. From vector field theory, (see 
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reference 3» section 1.7) » ve reedily establish in three-dimensional 
space. 


7^— — i-r— ■ -Ue6 (x - x’ ) 
It - Ct| 


(15) 


where 6 Is the Dirac delta function. (In two-dimensions, the multiplier 
in Fq. (15) becomes 2ir.) Furthermore, we know 

f^t(x)6(x - x')dT ■ f(x') (l6) 

where f is an arbitrary function. Substituting Ekis. (13) - (l6) into 
(12), we obtain 


Utt^(x) ** -^gT^^dx' • ftdo* (17) 


For our pertxurbation potmtial flow probloa, the first ri^t side 
term in Eq. (17) vanishes identically, due to Eq. (9). This is per- 
haps the single key feature of this approach in practical aerodynamics 
analysis, since only integrals on the closure 3R need be evaluated, 
and only surface velocities and pressures are of practical usefulness. 



(Were the potential Eq. (9) to he the inhomogeneous Poisson equation* then 
quadrature throu^out R would he required in Eq. (17) eucd the method 
loses considerahle computational appeal. ) !Ihe lead term in the surface 
integral can he evaluated hy a limiting operation* using the normal 
velocity boundary condition* Eq. (lO)* and the fact that -»■ 0 as 

|xj -*• * . Eq. (17) then becomes* for incompressible perturbation poten- 
tial flow, using Eq. (13) (See reference 1* page lU) 


2tr<Kx) - /«-0(x* )^*' 7 ” ' — • “do* ■ -n(x) • [U i - U] (l8) 


The last term in Eq. (18) represents the specified velocity efflux boun- 
dary condition statement. 

Using various potential flow singularity fvmctions (e.g. line source* 
vortex* do\iblet* see reference U)* the integrand in Eq. (Ic can be formed 
and eve0.uated. This yields a full* large order* linear algebraic equation 
sysxem to be solved for the strength distribution, o(x), of the elemen- 
tary singularities* as 


[A] { o) - { f} 


(19) 


These procedures can be extended into the high subsonic subcritical* 


and the linearized small disturbance purely s\q>ersonic flow regimes* 


U’lng a coordinate tranafoivation (see reference 5)« Hov including 
moKntuB considerations (developed in detail in the next section) » 
it is well knovn that the small disturbance linearized perturbation 
velocity potential flow equation becomes 


[1 - >di*. 


XX 


-f 

yy 


= 0 
ZZ 


( 20 ) 


The boundary c<xidition remains Eq. (lO), and is the 

A 2 2 

number associated with U i. Since 3 = [l-M ]isa 

00 oo 

(20) can be returned to the form of the Laplacian, Eq. 
independent variable coordinate transfoxmation 


reference tfach 
constant, Eq. 
(9 ) , by the 


X = X 

y = By (21) 

z = 3z 


Hence, incompressible .flow results can be extended to appropriate Mach 
nuniber ranges \ising Eq. (21) and solution to Eq. (9). For supersonic 
flows, B becomes coiqplex and physical velocity conponents are extracted 
from the real part of the coBg>\ztational solution. No transonic flow 
results are available using this approach. 



FORMULATION OF THE NONLINEAR PROBLEM 


In a V8QT, it is remarkable the extent to tihich the solution of the 
Laplanlan of the incui^ressible flow perturbation potential function 
has been successfully examined in potentie^. flow aerodynamics. We seek 
to evaluate the more comnrehensive tensor fom of the governing equation 
for irrotational potential flow, so as to exti -^t a potentially useful 
Green's function solution fom, valid throu^out a Mach number range, 
sx>ecifically including transonic flow. For the time-varying, isoenergetic 
flov of an inviscid perfect gas, conservation of mass, momentum and energy 
is expressed as the differential constraints. 

P»t ♦ “ 0 

‘’^i.t ^ PVi.J " P’i * ° 

2 

c - P,j * 0 


In the last equation, c is the local speed of sound, defined in terms of 
a stagnation condition eis 


2 

c 



Y_ - 1 


2 


U U 
k k 


(23) 
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The five equations (22) can he eoQvenien.t]y eohbined into a single 
equation, hy ei^anding the divergence in the first, substituting the 
density gradient by pressure fr«B the fifth, and eoabining into the 
reaiainder. In esqpanded fora, this yields 



- U U U 
i J i.J 


• I (U^).t * in P.t 


(2l*) 


In Sqs. (22) - (23), the subscript coana J denotes the gradient opera> 
tor on Xj , end the subscript conaa t denotes the partial derivative on 
time. For the additional constraint of irrotationality, lising Eq. (3), 
£q. (22t) takes the c<»ipact form 






+ In p,^ 


(25) 


The boundary condition for Eq. (25) remains Eq. (5). For steady flow, 

2 

the right side vanishes identically. Dividing through by c yields 
the familiar form 



-2 

c 






(26) 
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Ixi. (26) is the fundamntal fora for vbich ve seek to evalTiate 
existence of solution foras. We shall also have occasion to explore 
use of Eq.. (25 )» as a consequence of difficulties associated with 
Sq. (26). We cannot proceed directly to a Green’s function forauila- 
tion with £q. (26} since* a) it does not exist in divergence form, 
and/or, b) the coefficients of the Madi tensor, neither con- 

stants nor functioos of the spatial coordinates, but instead hif^ily 
non-linear functicsis of the dependent variaible. 

Eq. (26) explicitly displays the eonplete mixed elliptic-parabolic- 

hyperbolic character of potential flow. For exaiiq>le , for i = j , the 

2 

corresponding term in is of the fora [1 - m{i)]. where is the 

local Ma<^ nusber associated with the velocity component parallel to 
the i coordinate, i.e., Hence, for subsonic flows, is non- 

negative throughout R and Eq. (26) is elliptic. We therefore require 
knowledge of 4 or * n on the c<»q>lete closure 9R for a well-posed 
problem. For siperscmic floira, at least one of the diagonal entries 
in M. . is negative, and Eq. (26) displays a hyperbolic character. The 
well-posed solution requires knowledge of 4> and ^4 ^ n on a non- 
characteristic closure seffsent. In the intermediate range of primary 
interest (i.e., transonic), a diagonal coefficient in will be 
bounded arovmd zero, Eq. (26) will display a globally elliptic character 
with a priori unknown interior regions of hyperbolic and/or parabolic 
character. In this instance of major interest, we may seek to employ 
the transient description of Eq. (25), which can bring in an additional 
hyperbolic character to uniformly siqperimpose on the solution. 



Eqs. (26) and (25} are essentially Intractlble » tut sinileo'ly 
describe a much more general flov field than is of aerodynamic interest. 
Therefore, ve seek to eii^loy the simplifying small disturbance approxi** 
mation, but to retcdn the intrinsic tensor character while doing so. 
Hence, assume again that the external reference flow is parallel to 
the X axis (i.e., x^), and express the local velocity vector as a 
perturbation on of the form 


U. 




- <l>i 


(27) 


In Eq. (27), is the Kronecker delta which is zero except when k = 1. 
The reference state for the local sound speed can be efficiently 
modified to U^, yielding for Eq. (23) 


2 

c 




k k 




(28) 


As a first step towards tractibility, under the small disturbance approxi- 
mation for all subsonic and moderate supersonic (but explicitly not trans- 
onic) flows, Eq. (26) using Eqs. (27) - (28) will agedn yield Eq. (20). 
Rewritten in tensor form, Eq. (20) is formsHy identical to Eq. (26) 


with the definition 



(1 - 0 0 

0 10 

0 0 1 




(29) 


The last form empuaslzes that is only a scalar modification of the 
Kronecker delta in principal coordinates. With this identification, 
Eq.. (26) remains a linear partial differential equation but is no 
longer self-adjoint. To obtain the Green's function form, analogous 
to Eq. (12), we must identify the adjoint operator M operating on the 
Green's function g as (see reference 7» section 10 and reference 8, 
appendix I) 


Mg = (M^jg) 


(30) 


Denoting the differential operator on 0, Eq. (20), as L, form gli(fr, 
i.e. , 




( 31 ) 
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Inserting Eg. (31) into the divergence theorem, Eq. (ll), and noting 


that 


is constant, obtain 


* /3R®”l3*-l''3'“ - /9H<«“l3 > ’3'^l“ 

< 32 > 


Noting that the last term in Eq. (32) is £q. (30), and rearranging, ve 
obtain the tensor form of Green's third identity as 


- ♦MgldT (33) 


In writing Eq. (33), advantage was taken of rearrangement of dummy 
(repeated) subscripts. The first left term of Eq. (33) vanishes every- 
where within R, see Eq. (26). For boundary conditions. 









( 31 .) 


using Eq. (lO). Noting that H. . . vanishes* Eq. (33) can he rearranged 

1 J fi 

to the final form* analogovis to the scalar Green's theorem* Eq. (17)» 
as 






(35) 


Eq. (35) is the final form* valid for both subsonic and purely 
supersonic flows under the small disturbance approximation. A few 
additional comments are warranted* however. For the si^ersonic flow 
case* use of Eq. (35) rather than the Goethert luLe transformation of 
the incoBipressible flow solution* Eq. (9) mi^t yield improved solution 
accuracies since the boundeuy conditions are applied to the physical 
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aerodsmaalc shape (rather than the transfonned shape, see reference 6, 
pa^ 8). Caution Is required, houever, in eaplc^ng singularity methods 
to evaluate terms in Eq. (33)< For the supersonic linearized case 
imder consideration, Eq. (26) is hyperbolic, and g is more properly 
called a characteristic function (see reference 7, page 53). 

A hyperbolic equation does not admit to Isolated singularities, 
but instead every singularity is continued along a characteristic. 

Hence, theoretical evaluation of the character of g, end the all-in^rtant 
zones of silence ahead of a Mach cone, is required to assure proper 
application of Eq. (35) for > 1. No such problems are encountered 
for the high subsonic subcriticed case, and g takes on the familiar 
character of the Green’s function. We defer this presentation pending 
establishment of the transonic flow equation. 

The development leading to Eq. (35) required simplification of the 
basic equation using the small disturbance supersonic approximation. 

These linearizing assumptions are mutually exclusive to the transonic 
flow region of interest, and it is necessary to develop an alternative 
approach. As the first step, we seek to establish Eq. (26) in diver- 
gence form so as to employ the divergence theorem. Hence, using the 


chain rule 



I I - - I 1 


I 


Multiplying M. . ty c for convenience, the last term in E<i* (36) is 

* j 




i 




(37) 


Inserting Eq. (28) for and Eq. (27) for into Eq. (37)» we obtain 
equivalently 


'”*ij ■ 


ij 




" ■ *’i' 2 '*13 ■ ’^* 11*31 * 

c 


'S 

Im. W 


4>.i 

- *ii*^ji 


. Y..- 1 


♦’ll ♦’i^’n' 

H -1-^ - 2«- ^k3 - -H> 

L c 00 j ^ > 

m OO « 


(38) 


The last form in Eq. (38) utilizes the definition of 
Eq. (28), we have 
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( 39 ) 


r2\ 

2 

J 


-1 



We nov vish to evaluate Eq,s. (38) - ( 39 ) under the small dlstur- 
hance supersonic approximation, vhereupon ve assume 




(ho) 


Since ^ 1 in the transonic and lov supersonic range, this implies 

2 . Substituting for c^ in Eqs. (33) - (39), we observe that 

2 2 —1 

(c»/e ) ~ 1, and that all the non-constant terms within the outer 

00 

square bracket in Eq. (36) vanish Identically. Since the derivative 
of vanishes, we have established the conditions 

appropriate for Eq. (38) to vanish entirely. Hence, under the assump- 
tion of Eq. (UO), which restricts validity to a region around » 1, the 
fundamental potential flow equation, Eq. (26), can be eqxiivalently 
written in divergence form as 




(Ul) 
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This Is the desired result. It is of crucial iDportance, for 
practical cos^utational usefulness of the Green's function method, that 
the goreming partial differential equation he homogeneous, so as to 
avoid integrations over the entiz« solution domain R, see Eq. (17). 

Note that remains the complete tensor including off-diagonal 
entries, and thus encmi^asseB the linearised ca^^resslhle form, Eq. (29). 

Using Eq. (hi), and identifying G(x^,x^) as the nev tensor Green's 
function to he established, ve can obtain the desired solution form 
using the divergence theorem and Eq. (12) as 




The last form in Eq. (U2) is achieved recognizing that ”ij symmetric. 

The first left side term vanishes identically throughout R due to Eq. (hi). 
From the second left side term, ve vlsh to extract the field point value 
of potential, i.e., 9(x^). The form is again given hy Eq. (13); ve seek 
to transform the integral as given hy Eq. (l6). The desired form of the 
Green's functicm is (ref. 9) 


G(x^. xj) - I A 


m 


vhere. 


I M^^(x, - x')(x - xl) 
. . J»1 1 1 J J 


ikk) 


Is the inverse (matrix) of the Mach tensor • ^4* (26), and 
T^A ■ * 0. £q. (M») presents the necessary generalization of the 

geodesic distance r, Eq. (lU), associated with the scalar formulation. 
For the term of interest in Eq. (1*2), we obtain using Eqs. (U3) - (1»U) 






(Il5) 


Since , to the order of approximation of this development , derivatives 
on vanish, Eqs. (38) - (UO), so will they vanish on Performing 

the indicated differentiation on x. and x, in Eq. (1^5) will bring 



into the nyrorator while leaving [Vr^] in the den^^nator. Can«* 
cellation yi-rids tfc ' coefficient t the Krwtecker delta. 

The indices of differentiation then become contracted over i and J 
yielding the familiar evaluation given in Eq. (l$). Using Eq. (l6), 
ve then obtain 


/g«(xj) . -Uxt(Xj) 


(W) 


The Green's function solution fora for the non-lineeu transonic flov 
potential equation then becomes 


with the Green's function, given by Eqs. (U3) - (UU), 

Equation (1»7) is the terminal fora ve seek for transonic flov. 

As before, the boundary condition statement, Eq. (3) or (10), on 
velocity efflux over 6B win be used to eval\iate the first surface integral 
t -^ra. For transonic flow, it is highly iiq>robable that will introduce 
an additional singularity into the integral, since non>zero velocity 
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boundary ccmditions are specified veil avay from the transonic region. 
The seccmd integrand in Eq. (U7) is again en^loyed to fill out the 
influence coefficient matrix. The singularity that arises as x^^-*- x^ 
(where the superscript bar constraios x^ to 5 r) may be isolated by 
encircling x' with a small hemispherical surface of radius e. Pro- 
ceeding to the limit as e vanishes » a useful foimula is 


lim 



x! )n,da' 
1 J 


•1 

2 


(U8) 


Singiilarity distribution functions that are not constant over panels of 
the discretization, Eq. (1*7), will require individual evaluation of 
limiting behavior of the form of Eq. (1*8). 

While Sq. (U7) represents the Green's function solution form for 
transonic flow in completeness, several najor application-oriented 
features remain to be resolved that will undour,edly req\iire some 
computational evaluation. First, when the flow goes locally sonic, 
will vanish identically. If were a diagonal tensor, the determi- 
nant woxxld also vanish, and the inverse tensor required for the 
Green's function would becraue singulaur. However, in the present 
formulation, M^j is not diagonal but possesses the off diagonal ele- 
ments , In particular, miist be sufficiently large at 

the sonic point to keep det|M. ,| non-zero for the algorithm to remain 



stai>le. Of course, this is a function of the particular ge(^try, 
and probably requires a numerical evaluation. Downstream of this 
critical location for some distance, the flow will be locally super- 
sonic on the aerodynamic surface. In this region, will be negative 
and £q. (^7) will display the mixed hyperbolic character for which the 
tensor Green's function formulation is unproven, as discussed for Eq. 

(3^). Not all source points can exchange information with a field 
point in this region, as a manifestation of the zone of silence. Analyses 
of the transonic problem in this region using finite difference (ref. 10) 
and finite elenent (ref. 11 ) solutions to the field problem description 
have altered their difference kernel, and discarded negative "diffusion” 
coefficients, respectively to accurately model the locally siq>er8onic 
flow regions. Whether this can be suicoapllshed using a pure domain- 
closure procedure remains to be verified. Finally, for moderately high 
sTibsonic Incident flow Mach numbers, the s\q>ersonlc flow region will 
typically return to subsonic flow through a shock. The shock can be 
assumed to intersect the surface normally; since the flow parallel to 
each side of the shodc remains unlfom, the normal efflux boundary con- 
dition statement, Eq. (5), remains unaltered. However, so&e means will 
be required to accurately register the location of the shock, and to 
therein employ the Rankine-Hugoniot relations rather than Eq. (h7), 

(see reference 12, page 13). For solution of the field form of the 
transonic flow, artificial viscosity is typically employed to smear the 
shock over several mesh points and thereby eliminate the discontintiity. 


An alternative to the presented fonaulation exists which nic^t 
alleviate these critical flow recognition prohlens. Since the coefficients 
of the Mach tensor are non-linear functions of the potential function, 
an iteration procedure will be required to establish solution of Eq. (UT), 
in a manner corresponding to [A] becoming [A(o)], Eq. ( 19 ). The al- 
ternative would be to reformulate the description in terms of the trans- 
ient flow, Eq. (25). The addition of time as an independent vetriable 
can cast a uniform parabolic or hyperbolic character onto the solution. 
Using a four- vector approach (see reference 3 , chapter 6 ), a time-retarded 
Green's function could possibly be developed to account for a finite 
speed for signal propagation. For a scalar differential equation, the 
retarded Green's function form is 


g(x^, t; x’ , t' ) = 


|Xi - 


X. •“ X 

6{t' + — - t) + A 


ih9) 


where c 
at (x^. 


is the scalar signal propagation speed, the scalar equivalent of 
The causal description in Eq, (U9) indicates a signal is recorded 
t) due to a distiarbance originating earlier at (x^, t') where 


t' 


= t 




c 


(50) 
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For a peurabollc fom of the texiaor eq\iation, Eq. (25), a tf.me-retarded 
Green’s function is given as (see reference 9) 


G(xj^, t;x^, t’) 


3/2 


(t - f) 


[-R^A(t > f)] 


(51) 


2 

where R remains the geodesic distance, Eq. (UU). The equivalent ten- 
sor form for a hyperbolic equation is unknown. Its usefulness appears 
limited however, since in four vector form (see reference 13 for the 
completely linearized form), Eq. (25) is a generalized Poisson equation 
The nonhomogeneity is (in p),. ; on several occasions we have noted 
the practical usefulness requirement that the parent partial diffe- 
rential equation description be homogeneotis . In the final approach to 
steady-state, this term of course will vanish. However, should its 
retention be required to get a transient solution started, then in 
Green’s third identity, Eq. (U2), we would have the fotir-vector equiva- 
lent 




(52) 


Therefore, the formvilation requires evaluation of an integral over the 
entire domain R, and the key feature of only svirface paneling require- 
ments is lost. From a user standpoint, as well as cconputer pz-ogram 
design, this appeeurs totally unacceptable. 


CONCLUSIONS 


The existence of Green's function solution formulations for non- 
linear potential flov has been determined » for the tensor form of the 
governing differential equation system. For the smadl disturbance, 
completely-linearized approximation, and for subsonic and purely super- 
sonic flow, the Green's fimction solution form is given in Eq. (35)* 

For the transonic Mach number range, the general tensor equation is 
cast into divergence form, Eq. (Ul), using the small distuirbance approxi- 
mation. Subsequently, the Green's function is established for this ex- 
plicitly non-linear tensor differential equation, Eqs. (Us) - (Ul|), 
and the Green's function solution form established for transonic flow 
predictions, Eq. (U7). Summary ccmments are made regarding anticipated 
problems with application of Eq. (U7) to critical flow configurations. 

A plausible edternative solution procedure involving time-retarded 
Green's fxmctions is briefly described, although practical programming 
considerations contraindicate its utility. 
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